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In this paper the the effect of Dzyaloshinskii-Moriya interaction and anisotropy on the Entangle- 
ment of Heisenberg model has been studied. While the anisotropy suppress the entanglement due 
to favoring of the alignment of spins, the DM interaction restore the spoiled entanglement via the 
introduction of the quantum fluctuations. Thermodynamic limit of the model and emerging of non- 
analytic behavior of the entanglement have also been probed. The singularities of the entanglement 
correspond to the critical boundary separating different phases of the model. The effect of gapped 
and gapless phases of the model on the features of the entanglement has also been discussed. 
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I. INTRODUCTION 

In recent years, the notion of entanglement has re- 
ceived much attention in quantum information theory 
due to its notable features in developing the idea of quan- 
tum computers and other quantum information devices. 
Entanglement is a purely quantum correlation without 
classical counterpart 1 - and has been realized as a crucial 
resource to process and send information in novel ways 
such as quantum teleportation, supercoding and algo- 
rithms for quantum computations 2 -. Concerning the cor- 
relation content of the entanglement, states of systems in 
condensed matter physics may deserve the investigation 
of entanglement as a unique measure of quantum corre- 
lations. The interest will be intensified when we consider 
the relation between entanglement and quantum phase 
transition where a drastic change in the ground state 
of the system occur—. This change will occur at zero 
temperature where all thermal fluctuations get frozen 
and only surviving quantum fluctuations drive the phase 
transition. In the past few years the subject of many 
activities were to investigate the role of entanglement in 
the vicinity of quantum critical point for different spin 
modela 4 ! 5 ! 6 ! 7 - 8 - 9 - 10 ! 11 ' 12 - 13 . 

Spin models provide not only a test ground for above 
issue but also as a play-ground for implementation of may 
quantum information protocol o i 15 . Among them are 
Ising model in transverse field (ITF) and XXZ models. 
Despite of their simple Hamiltonian, low energy behav- 
ior of many systems can be captured through them. Ising 
model in transverse field has the benefit of exact solvabil- 
ity through mapping to the free fermion model^. Such 
solvability provide the possibility to test the behavior of 
entanglement and its scaling close to the critical point of 
the system and perform a finite size scaling as has been 
done in the seminal work of Osterloh and his coworkers 4 . 
The scenario is different in the XXZ model where the 
entanglement between the two nearest-neighbor sites de- 
velops a maximum at the isotropic point (A = 1) with- 
out any singularity in its first derivative 1 ^ which vanishes 



at the critical point A = 1. However, the block-block 
entanglement 1 ^ of the spin- 1/2 XX model with three-spin 
and uniform long range interaction shows a logarithmic 
and algebraic dependence on the size of block for differ- 
ent phases. Logarithmic divergences of the entanglement 
entropy is a general feature of all one-dimension critical 
systems where the coefficient of the logarithm is just the 
central charge of the underlying critical theory^. 

The scaling of entanglement close to the phase tran- 
sition and its connection to the universality class of the 
model can be further investigated through employing the 
renormalization group. This method as we will see in 
the next sections provides a rather analytic framework 
for treating the phases of the model even for those that 
are beyond the exact solution. In this stream the scaling 
of the entanglement govern by the critical exponent of 
the mode l 19 i 20 . However, The renormalization of quan- 
tum states has also been introduced in terms of matrix 
product states 21 . 

Both Ising and XXZ models can be supplemented with 
a magnetic term, the so called Dzyaloshinskii-Moriya 
(DM) interaction, arising from the spin-orbit coupling. 
Based on the symmetry aspects^ 2 -, It can be derived mi- 
croscopically as a linear correction to the standard su- 
perexchange mechanism™. The interaction has the form 
J2<ij> Dij.(Si x Sj) where the sum is over the pairs of 
spins. Some quantum antiferromagnetic (AF) systems 
are expected to be described by DM interaction with the 
underlying helical magnetic structures. Ising model with 
DM interaction was extensively studied^ 4 . DM interac- 
tion drive the quantum fluctuations resulting in a phase 
transition in the model. Critical point separates the an- 
tiferromagnetic and chiral phases. Derivative of the en- 
tanglement diverge at the critical point with the critical 
exponent of the model. 

In this paper we address the behavior of the entangle- 
ment in the XXZ model with DM interaction. Including 
the DM interaction makes the phase diagram rich with 
a critical line instead a single points. First we employ 
the quantum renormalization group to have a tractable 
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problem. Afterwards, the entanglement between degrees 
of freedom is treated through the renormalization. We 
will see that derivative of the entanglement becomes sin- 
gular at the phase boundary and its scaling correspond 
to the gapless and gapped phases of the model. The or- 
ganization of the paper is as follows. In the next section 
we briefly introduce the renormalization group approach. 
In the Sec (|HI[) we exemplifies the effect of anisotropy and 
DM interaction. Then in Secl|IV[) we turn on to discuss 
the scaling of the entanglement, and the Sec((V]) is de- 
voted to the conclusions. 



II. QUANTUM RENORMALIZATION GROUP 

The Quantum renormalization group presents a 
tractable version of treating quantum systems at zero 
temperatures. Through the renormalization the original 
model Hamiltonian is replaced by an effective one in cost 
of renormalizing coupling constants. In this way the orig- 
inal Hilbcrt space truncated to an reduced Hilbert space 
including the effective degrees of freedom. Getting ride 
of extra degrees of freedom gives rises to the flow of the 
coupling constants in the parameter space of the model. 
The version we employ to kill the degrees of freedom 
is Kadanoff's block approach since it is well suited to 
perform analytical calculations in the lattice models and 
they are conceptually easy to be extended to the higher 
dimension o 26 i 27 i 28 i 29 . In the Kadanoff's method, the lat- 
tice is divided into blocks in which the Hamiltonian is 
exactly diagonalized. By selecting a number of low-lying 
eigenstates of the blocks the full Hamiltonian is projected 
onto these eigenstates giving the effective (renormalized) 
Hamiltonian. 

The Hamiltonian of XXZ model with DM interaction 
in the z direction on a periodic chain of N sites is 



H(J, A) = ~ £ + °>l + x + A °i 



N 



+D(o?o» 1 -&»of ¥1 ) Z?,J,A>0 (1) 



We divide the chain into blocks each containing three 
sites. The block Hamiltonian (Hb = hf) of the three 
sites and its eigenstates and eigenvalues are given in Ap- 
pendix A of reference^. However, we give only the de- 
generate ground states since we need for them for eval- 
uation of entanglement and subsequent discussions, as 
follows: 



lift)) = 



f 



-(1 - iD)(A + q)\ HI) - 2[2iD + (D 2 - 1)]| TU>), (2) 



2(D 2 + 1)\ ITT) 



^/2 q ( q + A)(l + D 2 ) 
-(1 - iD)(A + q)\ UT) - + (D 2 - 1)]| TU») , (3) 




FIG. 1: Concurrence between first and third sites of a three- 
site model in terms of anisotropy for different values of DM 
interactions. 



where | f) and | j) are eigenstates of the a z Pauli operator 
and q = \/ A 2 + 8(1 + D 2 ). The projection operator of 
the ground state subspace defined by (Po = I 1T) (ipo | + 1 4 
)(V>o|), Where \ipo) and \ip' ) are the doubly degenerate 
ground states, | f|*) and | JJ.) are the renamed base kets 
in the effective Hilbert space. We have kept two states 
(\ipo) and IV'o)) f° r each block to define the effective (new) 
site. Thus, the effective site can be considered as a spin 
i. The effective Hamiltonian is similar to the initial one, 



Heff = ^E[<<+i + ^+i + A '^ + i 

i 

+D'(a?*V +1 -a?af +1 )] (4) 

where J' and D' are the renormalized coupling constants. 
The renormalized coupling constants are functions of the 
original ones which are given by the following equations. 



JV 



J' = J(-) 2 (l + D 2 ), 

q 



A' = 



l + D 2 



A 



-) 2 , D' = D.(5) 



The above RG equations show that there is a phase 
boundary A c = \/l 
p hase, A < \/l + D 2 

vTTc 2 — • 



D 2 that separates the spin fluid 
from the Neel phase, 



A > 



III. ENTANGLEMENT ANALYSIS 

Many measures of entanglement have been introduced 
and analyze d 30 i 31 ' 32 i 33 , but the most relevant to this work 
is the " entanglement of formation" . For a reduced den- 
sity matrix py of two qubits that arises after integrating 
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out other degrees of freedom, the entanglement between 
two qubits is evaluated as E = + — C 2 ), where 
h is a binary entropy function h(x) = —xlog 2 x — (1 — 
x) log 2 (l — x) and C denotes the concurrence^! defined 

as 

C = Max{Xi - A 2 - A 3 - A 4 , 0}, (6) 

where Afc(fc = 1, 2, 3, 4) are the square roots of the eigen- 
values in descending order of the operator Rij : 

Rtj = PijPij , Pij = {erf <g) oDpIj (of ® of). 

In this section we consider only a three site block and 
study the effect of DM interaction and anisotropy pa- 
rameter, i.e D and A, respectively, on the entanglement 
between two spins located on the sides of the block. To 
this end, let |^o) be the ground state of the block. By 
tracing the density matrix p = j^o)!^) with respect to 
the middle site of the block, the obtained reduced density 
matrix and Eq.® give the expression for the concurrence 
in terms of couplings D and A. 

For different values of DM interaction and anisotropy 
parameters, the plots of concurrence between first and 
third sites of the block C13 have been depicted in Fig.fl]). 
Consider first the case of D = 0. In this case the model 
becomes the known XX Z model. Large value of A im- 
plies the Neel state. Naturally this state is an product 
state without any entanglement between its constituents. 
As the anisotropy parameter reduces the quantum fluctu- 
ations arising from the transverse interactions have dom- 
inant effect and destroy Neel state. Indeed, the in-planar 
interactions drive the quantum correlations, i.e the qubits 
in the presence of quantum fluctuations are quantum cor- 
related. The main message of the Fig.([T]) is that the sup- 
pression of the entanglement can be compensated by tun- 
ing the DM interaction. In the absence of the anisotropy, 
the entanglement is insensitive to the DM value since 
both transverse interaction and DM term stimulate the 
quantum fluctuations. For nonzero value of anisotropy 
is clearly seen that the turning on the DM interaction 
restores the spoiled entanglement. The emerging of Neel 
state at large value of anisotropy and dominant quan- 
tum fluctuations at small value tempt to conclude that 
in the thermodynamic limit of the model there may oc- 
cur quantum phase transition with the critical boundary 
depend on the competition between the parameters in 
the Hamiltonian. We will address this issue in the next 
section. 



IV. THERMODYNAMIC LIMIT AND 
NON-ANALYTIC BEHAVIOR OF 
ENTANGLEMENT 

In this section we would like to see that how the quan- 
tum phase transition in the model, which can be signaled 
as unstable fixed point of RG equations, can be realized 
by examining the behavior of the entanglement. Indeed, 



A= 1 /T 
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FIG. 2: (color online) Representation of the evolution of en- 
tanglement entropy in terms of RG iterations at fixed value 
of anisotropy A = y/2. Two different behavior of the entan- 
glement at the large steps of RG correspond to the emerging 
phases of the model through the phase transition. 



the non-analytic behavior in some physical quantity is a 
feature of second-order quantum phase transition. It is 
also accompanied by a scaling behavior since the cor- 
relation length diverges and there is no characteristic 
length scale in the system at the critical point. As we 
already pointed out the renormalization group allows us 
to capture the thermodynamic properties of the model 
by considering a block of a few sites that is analytically 
tractable. In fact the global properties of the model enter 
a few sites through the renormalizing of coupling con- 
stants. We exploit this advantage to study the scaling 
of the entanglement in the model. Notice that in the 
n — th step of RG a system with size n^ +1 (ns is the 
number of sites in each block) describes effectively by 
a model consisting of only Jig sites with the renormal- 
ized of coupling constants. The case of the XXZ model 
has been extensively studied 2 *?, where the critical point 
A = f separates spin-fluid and Neel phases. However, 
for the present model the contribution of the planar DM 
interaction tune the critical point of the model due to 
involving the quantum fluctuations. Since the DM in- 
teraction doesn't flow, as is clear from the RG equations 
in Eq.©, it can be treated as a fixed parameter. Now 
we put the next step forward to see the evolution of the 
entanglement as the size of the system becomes large 
through the RG steps. 

Zero-step RG represents a three-site model that its en- 
tanglement studied in the preceding section. However, 
the first-step RG stands for a nine-site model which ef- 
fectively describes by a three-site model in the cost of 
renormalized coupling constants. In that case the entan- 
glement measures the correlation between effective de- 
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FIG. 3: (color online) First derivative of entanglement en- 
tropy and its manifestation towards divergence as the number 
of RG iterations increases (Fig[2]). 



FIG. 4: (color online) The scaling behavior of D m in terms 
of system size (N) where D m is the position of minimum in 
FigH 



grees of freedom. In each step of RG we can see the 
variation of the entanglement in terms of anisotropy pa- 
rameter with the fixed value of DM interaction. All these 
data have been shown in Fig. ([2]). In this figure we have 
set the A = \/2. It reveals that as the thermodynamic 
limit is touched via the increasing of the RG steps, the 
entanglement develops two rather different features. In- 
deed, there is a value for D = 1 that separates the dif- 
ferent features. This value is exactly the critical point 
of the model which is consistent with A c = yl + D 2 if 
we set A = \/2. Different features of the entanglement 
correspond to the emerging phases on both sides of the 
critical point. For anisotropy parameter larger than the 
critical point the Neel ordering dominates the phase of 
the model while for anisotropy parameter less than the 
critical value the increasing of the planar quantum fluc- 
tuations spoil any magnetic ordering. This feature is not 
a specific character of the model arising at A = ^/2. 
In fact for any value anisotropy A > 1 such behavior 
emerges with the only difference that the critical point is 
tuned into a new one. 

Further insight on the non-analytic behavior can be 
probed through the diverging of the first derivative of 
the entanglement at the critical point as long as the the 
thermodynamic limit is approached. Plots related to the 
derivative of the entanglement at different RG steps have 
been shown in the Fig. ([3]) . Each plot reveals a minimum 
which becomes singular as the critical point is touched. 
At the limit of large sizes of the model the singular behav- 
ior of the entanglement becomes more pronounced. One 
may wonder how such emerging singularity connects to 
the critical exponents or universality class of the model. 
To this purpose, we shall see that how position of min- 
imum D m i n and minimum value itself | \D mi „ scale 



with enlarging the size N of the system. Such a com- 
putation determines the scaling law of entanglement in 
one-dimensional spin systems and explicitly uncovers an 
accurate correspondence with the critical properties of 
the model. The position of the minimum (D m ) of 
tends towards the critical point like D m — D c — N~ 0A6 
which has been plotted in FigUJ Moreover, we have de- 
rived the scaling behavior of y = |gs|.D m versus N. This 
has been plotted in FigJS] which shows a linear behavior 
of ln(y) versus ln(N). The exponent for this behavior 
is | ~ N 0A6 . It should be emphasized that this 

exponent is directly related to the correlation length ex- 
ponent, v, close to the critical point. It has been shown 
in Ref.[l9j that | |g \ De ~ N 1 ^ and D m = D c - N' 1 ^. 

Singular behavior of 4^ corresponds to phase tran- 
sition for any value of DM interaction. It exhibits a 
singular curve at the transition points. The latter can 
be characterized by analyzing the derivative of entangle- 
ment for all values of DM interaction. As an example, in 
Fig. ([6]) the derivative of entanglement in a three dimen- 
sional view has been shown. Noticeably, the divergencies 
in the derivative are in perfect correspondence with the 
parameter's value at which the phase transition occur. 
The crack in the figure is just the critical line separating 
antiferromagnetic from the spin fluid phases. 

All above scaling functions hold for any value of 
anisotropy parameter as long as A > 1, which is a di- 
rect result of the fact that the parameter D doesn't flow. 
This means that the emerging DM interaction term in the 
model doesn't change the universality class of the model. 
Thus far, taken the derivative of entanglement with re- 
spect to anisotropy parameter, the singularity appears at 
the critical point. To get more insight about the role of 
DM interaction in the singularity of the entanglement, 
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ln(N) 



FIG. 5: (color online) The logarithm of the absolute value 
of minimum, ln(| dC/dA \ m in), versus the logarithm of chain 
size, ln(iV), which is linear and shows a scaling behavior. Each 
point corresponds to the minimum value of a single plot of 
FigS 
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FIG. 6: The crack appearing in the derivative of the entan- 
glement corresponds to the critical line of the model that sep- 
arates antiferromagnetic (right) and spin fluid (left) phases. 



it is convenient to plot the 4jy versus DM interaction 
as in Fig.©. Even at the high steps of RG no singu- 
larity detect. Observe that the pair (D = 1,A = \/2) 
stands for a point of singularity of derivative of entan- 
glement with respect to A as in Fig.©. However, there 
is no any signature of the divergence in latter quantity 
at this point when the derivative is taken with respect 
to D. This, again, verifies that the DM interaction have 
not any thing to do with the universality behavior of the 



FIG. 7: The derivative of entanglement ^ versus DM inter- 
action for fixed value of anisotropy A = \[2. Even in limit of 
high steps of RG no singularity is observed. 



model. 

Indeed, for A > vl + D 2 the long range behavior of 
the model fall into the the universality class of the Ising 
model that underlines the the appearance of the antifer- 
romagnetic long range order—. We emphasize that here, 
the critical line A c = \/l + D 2 from the Ising phase, 
i.e. D m — > D~. This directly comes from the fact 
that the Ising phase is a gapped phase. Approaching 
the critical point, the gap is closed as E g ~ \D — D c \" z , 
where z is the dynamical exponent. Since in the limit 
of large sizes of the system the critical point is touched 
as D c — D rn ~ AT -1 /", we are left with the result that 
the gap of the Ising phase in the proximity of the criti- 
cal point scales as E g ~ N~ z . Whenever A < VT+ D 2 , 
the model is gapless. This can be realized through a 
simple canonical transformation to the well known XXZ 
mode l 34 i 35 with the anisotropy A = ^^_ D -i ■ This implies 
that the model falls into a gapless spin fluid phase when 
A < 1. 

Through this paper we have only considered the the 
entanglement between two sites living on sides of a three- 
site block, i.e. the middle site has bees traced out through 
the reduced density matrix. We would like to emphasize 
that we could already consider the entanglement between 
first two sites of the block. The entanglement between 
first and second sites of the block has been shown that 
in Fig.©. At zero-step of RG that represent a three- 
site model, by increasing the DM interaction the entan- 
glement between sites reduce. This is in correspondence 
with behavior in Fig.© where the entanglement between 
first and third sites is increased by increasing the DM in- 
teraction. This is expected since when two parties get 
more entangled, they restricts their entanglement with 
third party and vis versa that is a reminiscence of the 
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FIG. 8: The entanglement between first and second sites of 
the block in terms of DM interaction at different RG steps. As 
before the anisotropy parameter has been fixed at A = \/2. 



monogamy property^ 6 , of entangled objects. For D < 1 
that corresponds to the Ising phase, the entanglement 
between two sites is shaved out. RG equations runs the 
anisotropy to infinity dictating the spins to align. In con- 
trary, for D > 1 that corresponds to the gapless phase, 
all plots behave independent of RG steps. Observe that 
for D < 1 neither first and second sites (C12) nor first 
and third ones (C13) in the large RG steps are entangled. 
This may not be surprising as in this limit the model is 
characterized by a polarized state. The situation is differ- 
ent for the gapless phase where the quantum fluctuations 
dominate the system suppressing the alignment of spins. 

If we were to take the derivative of plots, again the 
singularity reveals itself at the critical point. Although 
C13 and C12 present different behavior, they share in 
exhibiting the critical behavior of the model. 

V. SUMMARY AND CONCLUSIONS 

Condensed matter systems have received impetus from 
the concepts developed in quantum information theory. 
The central to this is the entanglement being a unique 
measure of the quantum correlations. In this stream 
we studied the entanglement in a one-dimensional mag- 
netic system in which many physical properties of real- 
istic complex materials can be understood through it. 
This model is the well known XXZ model supplemented 
by a magnetic term arising from the spin-orbit coupling. 
The phase diagram of the model is determined by the 



anisotropy and Dzyaloshinskii-Moriya interaction (DM) 
parameters. In a simple model consisting of only three 
qubits, the increasing of anisotropy parameter favors the 
alignment of spins antiferromagnetically yielding a prod- 
uct ground state without entanglement. However, tuning 
the DM interaction tends to build an entangled state and 
restores the spoiled entanglement. This reviving of en- 
tanglement can be understood via the fact that the DM 
interaction contributes the strong planar quantum fluc- 
tuations that pose the alignment ordering. 

The thermodynamic limit of the model realized 
through the rcnormalization group approach. The renor- 
malization group not only allows us to derive the critical 
points as well as phase diagram of the model but via im- 
plementation of the RG steps we are able to keep track 
the variation of the entanglement as the size of the sys- 
tem becomes large. RG equations imply that the DM 
interaction tunes the critical point, i.e there is critical 
line instead a single critical point. However, the uni- 
versality class of the model is unaffected in the pres- 
ence of the DM interaction, which can be clearly seen 
from both RG equations and the scaling we obtained 
for the entanglement. The role of the DM interaction 
can be well understood by analyzing and comparing the 
derivative of entanglement with respect to DM inter- 
action and anisotropy parameter. In the former case, 
even in the large RG steps, no singularity is observed. 
This can also be justified by mapping the model into the 
well known XXZ model using a canonical transforma- 
tion. The derivative of the entanglement diverge at the 
all points of the critical line. The line singularity corre- 
sponds to the phase boundary separating the antiferro- 
magnetic from the spin fluid phases. The singularity ac- 
companies by some scaling functions with a emerging ex- 
ponent that is related to the correlation length exponent 
close to the phase transition. We also verified that being 
gapped or gapless is relevant to the crossing the phase 
transition. Via the enlarging the size of the system, the 
singularity becomes more pronounced and touch the crit- 
ical line from the gapped phase (antifcrromagnetic). This 
phenomenon results the gap of the phase scales with the 
size of the system governed by the dynamical exponent. 
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